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Abstract. Zero-range effective interactions are commonly used in nuclear physics to 
describe a many-body system in the mean-field framework. If they are employed in beyond- 
mean-field models, an artificial ultraviolet divergence is generated by the zero-range of the 
interaction. We analyze this problem in symmetric nuclear matter with the to - h Skyrme 
^s ' model. In this case, the second-order energy correction diverges linearly with the momen- 

qq , turn cutoff A. After that, we extend the work to the case of nuclear matter with the full 

Skyrme interaction. A strong divergence (~ A 5 ) related to the velocity-dependent terms 
of the interaction is obtained. Moreover, a global fit can be simultaneously performed 
for both symmetric and nuclear matter with different neutron-to-proton ratios. These re- 
sults pave the way for applications to finite nuclei in the framework of beyond mean-field 
theories. 

"o ' 

1 Introduction 

Contact interactions are reasonable approximations for the real finite-range forces that can be em- 
ployed in cases where the interaction range is much smaller than the typical length-scale for the 

lO ' inter-particle distance. The main advantage of using zero-range interactions when treating many -body 

systems is that the equations to handle are usually simplified in this case. This is why contact inter- 
actions are extensively employed in many-body physics. Two examples of commonly adopted contact 
interactions are the zero-range effective Skyrme forces which are used in nuclear physics [ 1.2 1 and 
the contact interactions with coupling strengths depending on the s-wave scattering length which are 

^^ ■ employed for dilute atomic gases (3). 

In phenomenological effective interactions, such as Skyrme forces, a set of parameters has to be ad- 
justed. These parameters are commonly fitted to reproduce several observables at the mean-field level 
(for instance, in the nuclear case, binding energies and radii of some selected nuclei and the properties 
of nuclear matter); this means that these interactions are constructed to be used for mean-field -based 
calculations. When one goes beyond the standard mean-field theories, it is not obvious at all that the 
same effective interactions are still suitable to be used. An additional problem arises when these in- 
teractions are contact forces: In several beyond mean-field theories, an artificial ultraviolet divergence 
actually appears if the range of the interaction is zero. An ultraviolet divergence already appears in the 
mean-field-based Hartree-Fock-Bogoliubov (HFB) O or Bogoliubov-de Gennes (BdG) |5| models 
with zero-range forces in the pairing channel. Examples of beyond-mean-field models where ultravio- 
let divergences appear owing to the use of a zero-range interaction are: Higher- Tamm-Dancoff models 
where multiparticle-multihole configurations are introduced [6|, second random-phase-approximation 
models where not only ph matrix elements are included in the matrix to diagonalize owing to the pres- 
ence of two particle-two hole configurations [7 1 and models that describe the particle-phonon coupling 
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between individual degrees of freedom and collective vibrations |8 9). 

In nuclear physics, this divergence is usually treated by introducing a momentum cutoff A 11011111 or 
by applying the the techniques of dimensional regularization [ 121 . 

To summarize, when phenomenological zero-range interactions are used, two issues have to be ad- 
dressed if beyond-mean-field theories are adopted: 

1 . Choosing a strategy to handle the ultraviolet divergence associated to the use of a contact force; 

2. Deciding at which level the fitting procedure of the parameters has to be done and which constraints 
have to be included. 

In this article, we apply to a simple case a strategy to account for points 1) and 2). This strategy 
is quite general and can be applied to different cases and domains. We explore the possibility to fit 
phenomenological interactions at some beyond-mean-field level (so that these interactions are suitable 
to be used in beyond-mean-field models) by including the numerical momentum cutoff among the 
parameters of the interaction. 

The article is organized as follows. In Sec. [2] we consider symmetric nuclear matter with a simple zero- 
range density-dependent interaction that corresponds to the so-called to - tj, Skyrme model in nuclear 
physics. The energy per particle at the second-order beyond the mean-field approximation is derived 
analytically. In section[3] we consider nuclear matter with the full Skyrme interaction. Conclusions are 
summarized in Sec.|4] 



2 Second-order equation of state with the t Q - 1 3 Skyrme model 

In this section, we will consider symmetric nuclear matter and a zero-range density-dependent inter- 
action written as 

7(r 1 ,r 2 ) = <5(r 1 -r 2 )(f +|p a ) (1) 

where to, h and a are parameters. We introduce the strength g of the interaction: g(p) - to + 'j;p a ■ This 
choice of the interaction corresponds to the so-called to - h model that is a simplification of the usual 
Skyrme model where the non-local terms (i.e velocity dependent terms) and the spin-orbit part are 
neglected. Also, in nuclear physics, the finite-range Gogny force has a zero-range density dependent 
term. 
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Fig. 1. Direct and exchange first-order (upper line) and second-order (lower line) contributions to the total energy. 



The equation of state with the simplified Skyrme interaction V{v\,Vz) at the mean-field level for 
symmetric nuclear matter is calculated by considering the first-order diagram (plus the exchange term) 
displayed in the first line of Figfl] 
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The second-order correction is given by the diagram (plus the exchange term) plotted in the second 
line of Fig [T] 



& C -x -x -x v 2 (q) 
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e k 2 -q 



(3) 



In Eq.[3] the energies e are expressed as e* = |^-, where m* is the fermion effective mass. It is equal 
to the fermion mass m in the case of the to - h Skyrme model due to the absence of non-local terms. 
v(q) represents the interaction in the momentum space. 

A simple power counting argument shows the presence of ultraviolet divergence. Whereas k\ and £2 are 
limited by the Fermi momentum kf, q can be arbitrarily large and for large values of this momentum 
the second-order correction behaves like 

AEo,fd\^f. 

In other words, the second-order correction in symmetric nuclear matter with the simplified Skyrme 
interaction diverges linearly with the momentum cutoff A 1101111 : 



^(A A - 00) = ^L (_i 1 + 2 In 2) X (p) + A x(p) 
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where kf is the Fermi momentum and A is the particle number. 

In fig. |2|a), E/A + AE/A is plotted for different values of the cutoff A (from 0.5 to 2 frrT 1 ) and 
compared with the SkP-mean-field curve (solid black line). The correction AE/A is also plotted for 
different values of the cutoff A in Fig. EJb). One observes that the second-order correction causes a 
shift of the equilibrium point to lower densities. 
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Fig. 2. (a): E/A + AE/A as a function of the density p (in units of far 3 for different values of the cutoff A (in 
units of firT 1 ) compared with the SkP-mean-field equation of state (solid black line), (b): Second-order correction 
AE/A for different values of A. 



Eq. factually coincides with the equation of state obtained with the SkP Skyrme parameterization 
|[T3l where no contribution coming from the velocity dependent terms t\ and ti appears. In this work, 
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the SkP-mean-field equation of state is used to make comparisons with the beyond-mean-field equa- 
tions and is employed as our reference to perform the fit of the parameters. 

For each value of the cutoff A, we have performed a least square fit to determine a new parameter set 
SkP^ for the beyond-mean-field corrected equation of state E/A + AEjA so that the SkP mean-field 
EoS is reproduced. The corresponding curves obtained with the adjusted parameters are shown in Fig. 
[3] for different values of A. In the inset of Fig. [3] the refitted EoS are plotted and compared with the 
SLy5 mean-field EoS (solid line). The refitted parameters are listed together with the saturation point 
inTableQ] 



Table 1. From the second line, columns 2, 3 and 4: parameter sets obtained in the fits associated with different 
values of the cutoff A compared with the original set SkP (first line). In the fifth column the x 1 IN- value (x 1 divided 
by the number of fitted points) associated to each fit is shown. In column 6 the saturation point is shown. 





t (MeV fm 3 ) 


t 3 (MeV fm 3+3a ) 


a 


X 2 /N 


Po (fm- 3 ) 


E/A(po) (MeV) 


SkP 


-2931.70 


18708.97 


1/6 




0.16 


-15.95 


A = 0.5 fnr 1 


-2352.900 


15379.861 


0.217 


0.00004 


0.16 


-15.96 


A = 1 firT 1 


-1155.580 


9435.246 


0.572 


0.00142 


0.17 


-16.11 


A = 1.5 fnr 1 


-754.131 


8278.251 


1.011 


0.00106 


0.17 


-16.09 


A = 2 fnr 1 


-632.653 


5324.848 


0.886 


0.00192 


0.16 


-15.82 




0.05 



0.15 
P (fm" 3 ) 



0.25 



Fig. 3. Refitted second-order-corrected equations of state compared with the reference SkP-mean-field equation 
of state. 



3 Second-order equation of state with the full Skyrme interaction 



As it was done in section|2] we treat the equation of state of nuclear matter by adding the second-order 
contribution to the first-order mean-field energy. In this section, we use a standard Skyrme interaction 
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such as SLy5 in its complete form: 

V(ri,r 2 ) = f (l + XoPo-Mri - r 2 ) + $-(1 + x I P tr )[k' 2 5(r I - r 2 ) + <5(ri - r 2 )k 2 ] + f 2 (l + x 2 P IT )k' • <5(n - r 2 )k 



+ |(1 + Xi P IT ) P a p-±I?. Wr, - r 2 ) + iWoicri + cr 2 ) ■ [k' x <5(n - r 2 



)k]. 



In Ref. 1 1 1 ], the energy per particle in nuclear matter calculated at the second-order has been derived 
analytically using the general Skyrme force of Eq. [5] Its asymptotic behaviour for large values of 
the momentum cutoff A shows that the divergence is strong (~ A 5 ) and that this divergence is much 
stronger than the linear divergence of the to - t 3 model in Eq. [4) 



AE (2) * , (k F 

— - — (5,p, A — > oo) = ai(5,p) A + a 2 (6,p) A + a 3 (6,p) A + a 4 (6,p) + O I — 



(6) 



The last term a^iS, p) in Eq.|6]has been calculated explicitly in Ref [ 12 1 for symmetric nuclear matter 
(6 = 0), asymmetric matter and for pure neutron matter (6 - 1). 

In the upper panel of Fig. |4] we plot, as an example, the second-order EoS of symmetric nuclear 
matter obtained for different values of the cutoff A (see legend), from 0.5 up to 2 fin -1 . The different 
equations of state are calculated by using the SLy5 Skyrme parameters and are compared with the 
reference mean-field SLy5 EoS (solid line in (a)). In (b) the second-order correction is plotted for the 
same values of the cutoff A. 
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Fig. 4. (a) Second-order EoS for different values of the cutoff A and (b) second-order correction for symmetric 
nuclear matter calculated with the SLy5 parameters. The SLy5 mean-field EoS is also plotted in (a) (solid line). 



(5) 



In Ref. ifTTI . the problem of ultraviolet divergence has been solved in the following way. For each 
value of the cutoff A, we have performed a chi-squared minimization to determine the new sets of 
parameters for the beyond-mean-field corrected equation of state E/A + AE/A, All the parameters are 
kept free in the adjustment procedure. The minimization has been performed on 15 equidistant points 
using the following definition for the^ 2 : 



X' 



N 



-K 



N 'E-E A 2 



AEj 



(7) 



The errors or adopted standard deviations, AEj, are chosen equal to 1 % of the reference SLy5 mean- 
field energies E- ure f. This choice is arbitrary since we are fitting a theoretical EoS where a standard 



EPJ Web of Conferences 

deviation for this quantity has not been estimated. However, the magnitude of the y 1 defined in Eq. ?? 

has a clear and reasonable meaning: if it is smaller or equal to one, the reference EoS is reproduced 

within one standard deviation, i.e., within a 1% average error by our second-order EoS. 

In this work, a unique and global fit has been done to readjust the three mean-field plus second-order 

EoS for symmetric, asymmetric (6 - 0.5) and pure neutron matter to reproduce the corresponding 

SLy5 mean-field curves. 

The obtained sets of parameters are presented in Table |2] In Fig. [5] the three refitted EoS are plotted 

as function of the density p for different values of the cutoff A. The value of the saturation density is 

equal to 0.16 firT 3 for the four values of A. 



Table 2. Parameter sets obtained in the fit of the EoS of symmetric, asymmetric (6 = 0.5) and pure neutron matter 
for different values of the cutoff A compared with the original set SLy5. The standard deviation, cr, estimated for 
the different parameters is also given. In the last columns the x 2 and the incompressibilty K m values are shown. 
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0"ti 


0~h 


o-„ 


a-*) 


°**i 


°~X 2 


°"*3 


Wa 








(MeV fm 3 ) 


(MeV fm 5 ) 


(MeV fm 5 ) 


(MeV fm 3+3 ") 














(Mev) 


SLy5 


-2484.88 


483.13 


-549.40 


13736.0 


0.778 


-0.328 


-1.0 


1.267 


0.16667 




229.90 


/\(fm- j ) 




















X' 




0.5 


-2022.142 


290.312 


1499.483 


12334.459 


0.481 


-5.390 


-1.304 


0.880 


0.259 


0.411 


236.36 




0.49 


0.212 


1.75 


4.5 


0.001173 


0.00657 


0.00020 


0.001632 


0.000280 






1.0 


-627.078 


83.786 


-971.384 


186.775 


3.428 


-1.252 


-1.620 


200.360 


0.338 


0.540 


230.52 




1.668 


0.2740 


0.782 


0.078 


0.00260 


0.01927 


0.00026 


0.082 


0.000314 






1.5 


-743.227 


112.246 


-42.816 


5269.849 


1.013 


3.478 


-2.114 


0.189 


0.814 


1.733 


236.28 




0.306 


0.685 


0.2972 


5.4 


0.01415 


0.01309 


0.00519 


0.045037 


0.000784 






2.0 


-718.397 


573.884 


-497.766 


6179.243 


0.391 


-0.393 


-0.574 


0.785 


1.051 


1.313 


222.76 




0.343 


0.251 


0.261 


8.33 


0.005876 


0.001850 


0.000597 


0.017475 


0.00104 







The^- 2 in the global fit is composed by the three contributions and its final value is divided by three 
in order to make our different results comparable to one another. In other words, 



X' 



W* 



-■0)+X 2 (6 = 0.5)+x 2 (6=l)]. 



Globally, as one can see from the^- 2 values, the fit is of good quality. Specifically, these values are still 
less than 1 up to A = 1 firT 1 . Values between 1 and 2 (to be judged by considering the adopted choice 
of the errors in the expression of x 1 ) are found for larger values of the cutoff meaning that the fit is 
still good. Finally, we have estimated the standard deviation of the fitted parameters Ifl4l . This anal- 
ysis allows one to asses how well the used reference data together with the adopted errors constraint 
the parameters of our model. In particular, the standard deviation associated to such parameters are 
displayed in Table [2] 



4 Conclusions 



We have introduced a general method to handle the ultraviolet divergence generated by the use of a 
contact force in the evaluation of the second-order energy correction beyond the mean-field level. As 
a first step, we have considered symmetric nuclear matter with a simple zero-range density-dependent 
interaction, corresponding to the so called t^—tj, Skyrme model in nuclear physics. A cutoff regulariza- 
tion has been proposed to treat the ultraviolet divergence which is linear with the momentum cutoff A. 
We then proposed a fitting procedure such that the equation of state (EoS) including the second-order 
correction matches rather well the one obtained with the original SkP mean-field force. After that, 
we have extended our work to the case of nuclear matter with the nuclear Skyrme interaction. In this 
case, the velocity dependent terms of the Skyrme interaction have been included and the second-order 
equation of state shows a strong divergence (~ A 5 ). A global fit is performed simultaneously for the 
EoS of symmetric, asymmetric (6 = 0.5) and pure neutron matter at the second-order so that their 
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Fig. 5. Refitted EoS (global fit) for symmetric (a), asymmetric 6 = 0.5 (b), and pure neutron matter (c). The 
reference SLy5 mean-field curves are also plotted in the three panels (solid lines). 



corresponding SLy5 mean-field curves are well reproduced. These adjusted interactions display rea- 
sonable properties for nuclear matter. This opens new perspectives for future applications of this kind 
of interaction in beyond-mean-field models to treat finite nuclei. 
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